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THREE  DIMENSIONAL  CLOUD  DYNAMICS: 
PRELIMINARY  APPLICATION  TO  A  HANE  ENVIRONMENT 


I.  INTRODUCTION 

The  Defense  Nuclear  Agency  is  sponsoring  a  broad-based  research 
program  on  the  dynamics  and  phenomenology  of  high  altitude  nuclear 
explosions  (HANEs).  An  important  goal  of  this  research  is  to  understand 
and  model  the  late  time  nuclear  environment;  specifically,  the  intense 
field-aligned  striations  which  are  produced.  These  striations  are  regions 
of  high  electron  density  and  are  of  military  significance  because  they  can 
adversely  impact  communication  and  surveillance  systems.  Of  particular 
interest  to  DNA  is  the  development  of  a  model  which  characterizes  the  power 
spectral  density  of  nuclear  striations:  the  outer  scale,  the  freezing 
scale,  the  inner  scale,  and  the  associated  power  law  exponents  between 
these  scale  sizes.  This  information  provides  a  description  of  the  late 
time  nuclear  environment  which  can  be  used  in  propagation  codes  relevant  to 
military  systems.  '  03 

It  is  generally  believed  that  the  outer  scale  is  related  to  the 
macroscopic  size  of  the  striations  (or  clumps  of  striations),  while  the 
inner  scale  is  related  to  dissipative  effects  (i.e.,  diffusion).  The 
freezing  scale  length  is  associated  with  plasma  dynamics  at  scale  sizes 
between  the  inner  and  outer  scale  length  and  is  associated  with  a  'break' 
in  the  power  law  exponent.  This  type  of  power  spectral  density  is  observed 
in  naturally  occurring  structured  ionospheric  environments  (e.g., 
equatorial  spread  F,  high  latitude  blobs).  Since  it  is  believed  that  the 
fundamental  structuring  processes  are  the  same  in  both  the  ambient  and  HANE 
environments,  then  it  is  expected  that  the  power  spectral  densities  will  be 
similar,  at  least  qualitatively. 

Recently,  DNA  researchers  have  focussed  on  developing  a  model  for  the 
freezing  scale  length  (k^  ^ )  which  can  be  implemented  in  conjunction  with 
the  DNA  code  SCENARIO.  One  freezing  model  that  has  been  t eveloped  is  the 
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JAYCOR  freezing  model  (Sperling  et  al.,  1988).  This  model  attributes  the 

freezing  scale  length  to  collisional  viscosity  and/or  magnetic  viscosity 

depending  on  the  ionospheric  parameters.  This  model  emphasizes  the 

dynamics  of  the  plasma  transverse  to  the  magnetic  field  and  incorporates 

parallel  effects  in  a  qualitative  manner.  In  contrast  to  this  model, 

recent  studies  at  NRL  (Drake  and  Huba,  1986;  1987;  Drake  et  al.,  1988)  have 

shown  that  parallel  dynamics  (i.e.,  3D  effects)  can  play  the  key  role  in 

the  'freezing'  of  plasma  clouds  (or  striations).  In  particular,  a 

stability  criterion  has  been  derived  which  indicates  that  plasma  clouds  (or 

striations)  are  stable  against  large-scale  structuring  when  r^  >  r  ;  here 

r  is  the  radius  of  the  cloud  transverse  to  B,  and  r  is  a  'critical' 
c  ~  cr 

radius  which  is  a  function  of  the  physical  parameters  of  the  system.  The 
purpose  of  this  report  is  to  present  quantitative  estimates  of  rcf  for 
parameters  typical  of  barium  clouds  and  nuclear  striations  based  upon  the 
NRL  3D  model  and  to  discuss  the  implications  of  these  results  regarding  the 
freezing  scale  length. 

For  the  case  of  barium  clouds,  we  find  the  critical  scale  size  for 

marginal  stability  is  r  <  100  m,  i.e.,  barium  clouds  with  r  >  r  are 

cr  ~  c  cr 

stable.  Since  the  scale  size  of  'frozen'  barium  clouds  is  ~  100's  m,  we 
conclude  that  our  theoretical  results  are  consistent  with  the  observations: 
3D  effects  can  provide  the  physical  mechanism  to  'freeze'  barium  clouds  at 
the  observed  scale  sizes.  Although  this  result  is  very  significant,  it 
does  not  explain  the  actual  scale  size  of  the  striations.  We  suggest  the 
following  hypothesis  regarding  the  structuring  of  barium  clouds.  During 
the  initial  phase  of  a  barium  cloud  expansion,  the  cloud  evolves 
dynamically.  The  cloud  diffuses  along  the  magnetic  field,  and  in  the 
presence  of  a  neutral  wind  or  ambient  electric  field  it  steepens  on  one 
side.  Eventually,  the  steepened  'backside'  of  the  cloud  can  become 
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unstable  to  the  E  x  B  instability  and  the  cloud  structures.  During  this 
phase,  the  waterbag  cloud  model  used  in  our  3D  theory  is  not  appropriate, 


and  hence,  our  theory  is  not  applicable.  This  is  not  to  say  that  3D 
effects  are  not  important  during  the  onset  of  instability,  but  only  that 
the  quantitative  analysis  in  Drake  et  al.  (1988)  is  not  valid.  Subsequent 
to  the  initial  break-up  of  the  cloud,  the  smaller  clouds  (or  fingers)  may 
evolve  into  waterbag-like  equilibria.  If  this  is  the  case,  then  we  can 
apply  the  NRL  3D  model  and  ask  the  question:  Are  these  smaller  clouds 
stable  or  unstable?  If  they  are  stable,  as  our  analysis  indicates,  then  3D 
effects  can  account  for  'freezing'.  As  to  the  distribution  of  frozen  scale 
sizes,  this  would  be  determined  by  the  initial  break-up  of  the  cloud. 

Applying  our  model  to  plasma  parameters  believed  typical  of  nuclear 
striations,  we  find  that  the  critical  scale  size  is  in  the  range  rcr  =  10's 
m  -  100's  km.  The  reason  for  this  rather  broad  range  of  values  is  due  to 
the  broad  parameter  regime  considered.  Thus,  unlike  the  situation  for 
barium  clouds,  it  is  not  clear  that  the  initial  set  of  nuclear  striations 
will  be  stable  to  further  structuring  because  of  3D  effects.  Nevertheless, 
early  time  structuring  (t  <  few  min)  will  probably  play  an  important  role 
in  setting  up  a  distribution  of  stiiation  scale  sizes  which  may  affect  late 
time  phenomenology,  and  hence  should  be  more  thoroughly  investigated. 


II.  CRITICAL  RADIUS  MODEL 

The  plasma  and  field  configuration  used  in  Drake  et  al.  (1988)  is 

described  as  follows.  The  magnetic  field  is  in  the  z-direction  (B  =  B  e^), 

the  neutral  wind  is  in  the  x-direction  (V  =  V  e  ),  and  a  waterbag  model 

is  assumed  for  the  density,  i.e.,  nn  =  n  +  n,  inside  the  cloud  and  n~  =  n. 

J  0  c  b  Ob 

outside  the  cloud  where  the  subscripts  c  and  b  denote  cloud  and  background, 
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respectively.  The  cloud  is  assumed  to  be  circular  in  the  plane  transverse 
to  B  and  is  an  arbitrary  ellipse  in  the  plane  containing  B. 

The  equilibrium  potential  for  this  configuration  consists  of  two 
components:  a  polarization  potential  caused  by  the  neutral  wind,  and  an 
ambipolar  potential  caused  by  electron  pressure  parallel  to  B.  The 
polarization  potential  causes  the  cloud  to  move  in  the  direction  of  the 
neutral  wind  (albeit  at  a  reduced  speed).  The  relative  ion-neutral  slip 
velocity  is  the  driving  mechanism  for  the  gradient  drift  instability  which 
can  cause  the  cloud  to  structure.  The  ambipolar  potential  causes  the  cloud 
to  spin  about  the  axis  aligned  with  the  magnetic  field.  The  rotation  rate 
varies  with  position  along  the  length  of  the  cloud  (i.e.,  along  B),  so  that 
there  is  a  'shear'  in  the  cloud  velocity  in  this  direction. 

A  stability  analysis  was  performed  based  upon  this  equilibrium  and  a 
stability  criterion  was  derived  [Drake  et  al.,  1988].  The  marginal 
stability  criterion  for  the  cloud  is  given  by 


where 


T  = 


r 

c 


r  = 


V 


i_ 

V 

n 


1_ 

r 

c 


(1) 


(2) 


where  4>q  is  determined  from 

Wo(*0’  rc’  M’  Zc)  =  0  (4> 
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These  expressions  are  (52a)  and  (52b)  in  Drake  et  al.  (1988).  Here,  r£  is 

1/2 

the  radius  of  the  cloud,  zp  =  ,  Lz  is  the  half-length  of 

the  cloud  along  B,  cn  .  are  the  parallel  and  perpendicular  conductivities, 

1/2 

respectively,  [  af f/ =  (QeQ^/\>e\K)  where  Ga  is  the  cyclotron  frequency 
of  species  a,  v  =  v  .  +  v  and  v  a  is  the  collision  frequency  between 

species  a  and  0] ,  M  =  n^n^,  is  the  azimuthal  angle  for  mode 

localization,  Vq  is  the  cloud  speed,  and  Mq  is  the  azimuthal  mode 
frequency.  Clouds  are  stable  for  T  >  Tc  and  are  unstable  for  T  <  T  . 

This  stability  criterion  is  based  upon  two  distinct  physical  effects: 
convection  and  shear.  First,  the  condition  Wq  =  0  requires  the  mode  to  be 
non-propagating;  this  results  from  a  balance  of  the  diamagnetic  propagation 
of  the  mode,  and  the  azimuthal  background  plasma  flow.  If  this  condition 
cannot  be  met,  then  the  cloud  is  stable  because  perturbations  are  rapidly 
convected  from  the  unstable  'backside'  of  the  cloud  to  the  stable  'front- 
side'  of  the  cloud.  This  stabilizing  mechanism  is  discussed  in  detail  in 
Drake  and  Huba  (1986).  The  second  stabilizing  mechanism  results  from  shear 
stabilization.  As  noted  earlier,  the  azimuthal  flows  generated  by  the 
ambipolar  potential  are  not  uniform  along  the  length  of  the  cloud.  This 
sheared  flow  causes  the  perturbations  to  have  kz  t  0  which  gives  rise  to 
the  damping  of  the  mode.  This  point  is  discussed  in  detail  in  Drake  et  al. 
(1988). 

The  stability  criterion  can  be  simplified  by  considering  the  limits  z^ 
>>  1,  zc  =  1,  and  z^  <<  1  (Drake  et  al.,  1988).  We  find  that  (3)  can  be 
written  as 

2zc2  cos<*>0/(M  +  2)  zc  »  1 

T  =  2z  2  cos<J>n/ (M  +  3)  z  =  1  (5) 

ecu  c 

(2/fl)z  2  cos4>A/(2/it  +  Mz  )  z  <<  1 

c  0  c  c 
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We  further  simplify  (5)  by  taking  4>q  ~  0°  for  <  1  and  <|>q  ~  45°  for  z > 
1.  This  approximation  is  based  upon  exact  solutions  to  (1)  -  (4),  and  is 
used  in  order  to  develop  a  simple  model  for  the  critical  cloud  radius.  We 
obtain  the  following  marginal  stability  condition 

42  z  2/(M+2)  z  »  1 

c  c 

r  =  rc  £  342  zc2/2(M+3)  z c  =  1  (6) 

(2/n)z  2/(2/ji  +  Mz  )  z  «  1 

c  c  c 


However,  we  note  that  our  goal  is  to  determine  a  critical  striation 
radius,  *cr,  based  upon  (6).  In  order  to  do  this  we  proceed  in  the 
following  manner.  We  define  as  follows 


so  that 


ro 


lTe  +  T 


eB 


i}  1_  1_  (Vi172 
V„  L  UJ 


(7) 


r  =  rnz  . 

0  c 


(8) 


Substituting  (8)  into  the  marginal  stability  condition  (6)  we  arrive  at  a 
new  marginal  stability  condition 

42  z  /(M+2)  z  »  1 

c  c 

r0  =  rc0  s  2  zc/7(M+3)  zc  =  1  (9) 

(2/n)z  /(2/n  +  Mz  )  z  «  1 
c  c  c 


The  advantage  of  using  (9)  instead  of  (6)  is  that  Tq  does  not  depend  on  the 
transverse  radius  of  the  cloud  (or  striation). 


Based  upon  (9)  we  determine  the  critical  radius,  rcr,  as  follows. 

First,  we  plot  I*cq  vs.  z for  a  specific  value  of  M.  This  is  shown  in  Fig. 

1  for  M  =  5.0;  note  that  I^q  increases  monotonically  with  increasing  z^. 

In  calculating  T  n  we  use  (9)  for  z  >  10,  z  =  1,  and  z  <0.1.  The 

cO  c  c  c 

values  of  T  n  between  z  =  0.1  and  1.0,  and  z  =  1.0  and  10.0  are  calcu- 
cU  c  c 

lated  assuming  linear  interpolation  between  these  points.  We  find  this  to 

be  a  reasonable  approximation  based  upon  comparisons  with  exact  solutions 

of  Tcq  using  the  formulas  in  Drake  et  al.  (1988).  Second,  we  calculate  the 

value  of  Tq  and  plot  it  as  a  straight  line.  For  illustrative  purposes  we 

assume  Tq  =  0.5  in  Fig.  1.  The  marginal  stability  condition  is  determined 

by  rQ  =  rcQ.  Clouds  are  stable  for  TQ  >  and  unstable  for  TQ  <  T^. 

This  means  that  clouds  are  stable  for  z  <  z  and  unstable  for  z  >  z 

c  cr  c  cr 

1  /o 

where  z  =  2.2  in  Fig.  1.  Noting  that  z  =  (L  /r  )(<t./ctm)  we  then 
cr  c  z  c  1  II 

define  the  critical  radius  to  be 


.  ^  fV|1/2 
:r  z  U,,J 


(10) 


cr 


Thus,  clouds  are  stable  for  r  >  r  and  are  unstable  for  r  <  r  The 

c  cr  c  cr 

explicit  dependence  of  r^  on  the  various  parameters  (e.g.,  L z,  V^,  o^/o^, 
etc.)  is  described  in  the  Appendix. 


III.  QUANTITATIVE  RESULTS 

We  now  present  a  series  of  figures  which  plot  r  as  a  function  of 
different  physical  parameters  representative  of  barium  clouds  and  nuclear 
striations. 

First,  we  consider  barium  clouds.  In  Fig.  2  we  plot  r^  (meters)  vs. 

a. Jo,  for  M  =  10,  T  =  T.  =  0.1  eV,  B  =  0.3  G.  V  =  100  m/sec,  and  L  =5 
111  e  l  n  z 

km,  10  km,  15  km,  20  km,  and  40  km.  Barium  clouds  are  observed  to 
structure  when  ~  15  -  30  km  (Linson,  1972].  For  clouds  at  altitudes  ~ 
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A 

180  km  we  estimate  <r. . / cr .  ~  6.25  x  10  .  From  Fig.  2  we  see  that  r  <  10's 

II  1  6  cr  ~ 

m  for  Lz  <  40  km.  Thus,  barium  clouds  with  transverse  scale  sizes  greater 
than  10's  m  are  stable  to  further  structuring  by  the  gradient  drift 
instability  because  of  3D  effects.  This  is  consistent  with  observational 
evidence  which  suggests  that  barium  clouds  are  frozen  for  r£  ~  100's  m. 
However,  we  add  that  the  critical  radius  determined  by  our  theory  does  not 
predict  the  predominance  of  cloud  striations  with  r^  ~  100's  m.  Ve  discuss 
this  point  further  in  the  final  section. 

We  now  turn  our  attention  to  parameters  typical  of  nuclear  striations. 

3 

In  Fig.  3  we  plot  r  (km)  vs.  ct./ct.  for  L  =  200  km,  V  =10  m/sec,  T  = 
6  r  cr  II  1  z  n  e 

T^  =  0.3  eV,  B  =  0.3  G,  and  M  =  5.0,  10.0,  20.0,  and  100.0.  First,  we  note 

that  in  general,  r^  decreases  as  M  increases.  Second,  as  in  the  case  of 

barium  clouds,  there  is  a  relatively  strong  dependence  of  rcf  on 

For  cr./a.  >  4.0  x  10^  we  note  that  r  <  100's  m  for  the  values  of  M 

II  1  ~  cr 

chosen.  However,  for  more  collisional  striations  (i.e.,  a{\/ai  <  2.0  x 

10^),  the  critical  radius  for  low  to  moderate  M  values  is  >  10  km,  while 

for  high  M  values  is  ~  100's  m.  And  third,  note  that  r  becomes 

4 

independent  of  M  for  <  10  and  M  <  20;  for  these  parameters  we  also 

have  z  <<  1.  This  result  is  consistent  with  (10). 
cr 

3 

In  Fig.  4  we  plot  r  (km)  vs.  a. ./a.  for  M  =  10,  V  =10  m/sec,  T 
°  cr  II  1  n  e 

T.  =  0.3  eV,  B  =  0.3  G,  and  L  =  100  km,  200  km,  300  km,  and  400  km.  The 
l  ’  z 

values  of  r  are  similar  to  those  in  Fig.  3,  ranging  from  10  s  m  to  100's 

km.  We  note  that  r  decreases  as  L  decreases  as  expected.  For 

cr  z 

relatively  short  striations  (L  ~  100  km),  we  find  that  r  <  10  km  while 
J  z  cr  ~ 

for  longer  striations  (L^  ~  400  km)  that  r  <  100's  km. 

In  Fig.  5  we  plot  r  (km)  vs.  0../0.  for  M  =  10,  T  =  T.  =  0.3  eV,  B  = 
&  r  cr  II  1  e  1 

0.3  G,  L  200  km,  V  =  10^  m/sec,  6.7  x  10"  m/sec,  3.3  x  10^  m/sec,  and 
’  z  n 

2 

10  m/sec.  We  note  that  r  decreases  substantially  as  V  decreases  from 

cr  J  n 
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3  2 

10  m/sec  to  a  nominal  ionospheric  value  of  10  m/sec.  For  conductivity 

5  2 

ratios  >  10  we  find  that  r  <  300  m  when  Vn  =  10  m/sec;  this  value 

of  r  is  similar  to  the  barium  cloud  results, 
cr 

IV.  DISCUSSION 

We  have  developed  a  relatively  simple  procedure  to  determine  the 
critical  transverse  scale  size  of  an  ionospheric  plasma  cloud  or  striation 
based  upon  the  3D  model  of  Drake  et  al.  (1988).  The  critical  radius  is 
given  by  (9)  where  z is  determined  from  (8),  i.e.,  =  I^q.  We  find  the 

somewhat  surprising  result  that  'thin'  clouds  are  more  unstable  than  'fat' 
clouds  (all  other  parameters  being  equal).  That  is,  clouds  are  unstable 

for  r  <  r  ,  while  they  are  stable  for  r  >  r  .  Applying  this  model  to 

barium  clouds  released  at  ~  180  km  we  find  that  r^  <  100  m;  this  result  is 
consistent  with  observational  evidence  of  frozen  clouds  (or  fingers)  for 
~  100's  m  [Linson,  1972].  Applying  this  model  to  nuclear  striations  leads 
to  a  broad  range  of  rcj_  because  of  the  variability  of  the  physical 
parameters  of  the  system.  We  find  that  for  weakly  collisional  plumes 
(i.e.,  CT||/a^  >  10^)  that  r^  <  few  kms ,  while  for  strongly  collisional 
nlumes  (i.e.,  CT| ( / cr^  <  10^)  that  r  >  10  km. 

Although  these  results  are  interesting,  there  are  two  problems  which 
need  to  be  addressed  regarding  the  development  of  a  freezing  scale  model 
based  on  3D  dynamics.  The  first  problem  involves  the  observed  size  of 
frozen  barium  cloud  striations  (r  ~  100's  m).  This  scale  size  is  not 

described  by  our  3D  model.  An  explanation  for  this  is  the  following. 

During  the  early  stages  of  barium  cloud  evolution,  the  cloud  is  relatively 

diffuse  and  is  not  wel 1 -descr i bed  by  the  3D  waterbag  model  assumed  in  Drake 
et  al.  (1988).  The  cloud  is  not  in  equilibrium  and  evolves  dynamically 
with  its  'backside'  developing  a  steep  gradient.  The  'backside'  eventually 
becomes  unstable  and  the  cloud  'breaks-up'  into  smaller  clouds  (or 
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fingers).  We  speculate  that  the  observed  freezing  scale  size  is  determined 

during  the  initial  structuring  of  the  cloud.  The  newly  formed  smaller 

clouds  may  have  a  density  profile  similar  to  the  waterbag  equilibrium 

discussed  in  Drake  et  al.  (1988)  (or  may  evolve  into  such  a  profile),  and 

thus,  would  be  stable  to  further  structuring  because  of  3D  effects. 

Applying  this  reasoning  to  nuclear  striations  suggests  that  it  is 

crucial  to  understand  and  model  the  initial  onset  of  structure  in  the 

nuclear  environment  in  order  to  determine  the  physical  scale  sizes  of 

striations.  However,  unlike  the  barium  cloud  case,  it  is  not  clear  that 

the  initial  set  of  scale  sizes  will  be  stable  to  further  structuring 

because  of  3D  effects.  This  is  because  the  broad  range  of  physical 

parameters  (e.g.,  ai/a±>  ,  vn»  etc.)  in  a  nuclear  environment  leads  to  a 

broad  range  of  critical  scale  sizes  r  ;  as  shown  in  Section  III,  r  = 

cr  cr 

10's  m  -  100's  km.  It  is  likely  that  there  will  be  striations  formed  with 
rc  <  rcr  so  that  further  structuring  can  occur.  Nevertheless,  early  time 
structuring  (t  ~  few  min)  will  probably  play  an  important  role  in  setting 
up  a  distribution  of  striation  scale  sizes  which  may  persist  into  late 
time,  and  hence  should  be  more  thoroughly  investigated. 

For  example,  the  field-aligned  striations  which  formed  during  the 
Checkmate  event  had  transverse  scale  sizes  in  the  range  300  -  1800  m,  most 
of  them  were  in  the  range  500  -  700  m  (Chesnut,  private  communication). 

These  striations  were  produced  within  a  few  minutes  after  the  burst.  To 

O 

estimate  r  for  Checkmate,  we  consider  an  aluminium  plasma  with  n  =  10 
cr  1 

-3  -1 

cm  ,  Tg  =  Tj  =  .3  ev,  B  =  0.3  G,  and  =  1  sec  .  For  these  parameters 

we  find  that  a,,/a,  =  2  x  10^.  From  Fig.  3  we  note  that  r  =  2  km  (for  M  = 
111  °  cr  ' 

3 

100,  Vn  =  10  m/sec,  and  =  200  km)  which  is  larger  than  the  observed 

striation  radii.  However,  if  we  were  to  take  M  somewhat  larger  and  V 

&  n 
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somewhat  smaller  (e.g.,  M  =  200  and  V  =  400  m/sec)  we  find  that  r  =  300 

n  cr 

m.  For  these  parameters,  the  Checkmate  striations  would  be  stable  to 

further  structuring  by  the  E  x  B  instability;  the  important  scale  sizes  and 

relevant  power  spectra  would  then  be  determined  at  early  time. 

Up  to  this  point,  the  discussion  has  focussed  on  the  transverse  scale 

size  of  frozen  barium  or  nuclear  clouds.  Here,  frozen  refers  to  the 

property  that  the  plasma  cloud  is  stable  to  further  large-scale  structuring 

by  the  gradient  drift  instability.  This  leads  to  the  second  problem:  what 

is  the  relationship  between  the  frozen  transverse  scale  size  and  the 

freezing  scale  length?  The  freezing  scale  length  is  defined  as  the  scale 

size  at  which  there  is  a  transition  in  the  power  law  of  the  power  spectral 

-15  -25 

density  (for  example,  from  ~  k  to  ~  k  ).  The  details  of  the  power 

spectral  density  depend  not  only  upon  the  transverse  scale  of  the  cloud  or 

striation,  but  upon  the  distribution  of  transverse  scale  sizes  of 

striations,  as  well  as  the  density  profiles  themselves  (i.e.,  edge  effects) 

(Zabusky  et  al.,  1986;  Prettie,  1988).  Calculating  these  quantities  is  not 

within  the  scope  of  the  present  3D  model;  3D  numerical  simulations  will  be 

required  to  better  understand  many  of  these  issues. 

In  conclusion,  we  have  presented  quantitative  estimates  of  the 

critical  scale  size  (rcr)  of  ionospheric  plasma  striations  based  on  the  3D 

model  developed  at  NRL  (Drake  et  al.,  1988).  We  have  considered  parameters 

typical  of  both  barium  clouds  and  nuclear  striations.  The  relevance  of 

these  results  to  the  'freezing'  problem  is  described  as  follows.  If  the 

initial  set  of  striations  has  scale  sizes  such  that  r  >  r  then  they  will 

c  cr 

be  stable  to  further  large-scale  structuring  because  of  3D  effects,  i.e., 
they  will  be  'frozen'.  This  hypothesis  is  consistent  with  barium  cloud 
observations.  This  leads  to  the  important  conclusion  that  the  late  time 
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characteristics  of  striations  (e.g.,  scale  sizes,  power  spectra,  etc.)  are 
determined  at  early  time.  Thus,  a  thorough  investigation  of  early  time 
structuring  processes  is  needed.  Applying  this  idea  to  the  nuclear 
environment  suggests  that  the  structuring  of  the  debris-air  shell  (being 
studied  experimentally  by  the  NRL  laser  group),  as  well  as  the  possible 
structuring  of  the  deposition  region,  should  be  vigorously  pursued. 
Finally,  we  add  that  our  results  are  preliminary  in  the  sense  that  they  are 
based  upon  an  idealized  cloud  model.  To  better  understand  and  quantify  the 
nature  of  3D  effects  for  realistic  clouds  requires  3D  simulation  studies; 
such  studies  are  currently  being  conducted  at  NRL. 
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APPENDIX 


We  present  an  explicit  expression  for  rcr  as  a  function  of  the 

parameters  (e.g.,  M,  Vn,  Tg,  T.,  a\\/ai’  Lz’  etc->-  Based  upon  the 

1/2 

definition  z c  =  (Lz/rc) ( a^/ a^)  ,  we  rewrite  T ^  (Eq.  (9))  as  follows 


L  ,a,,  1/2 


U  /W  I  V  A  / 

J2  f  MO  («  ♦  2)' 


-c  "'ll 


r  «  L 
c  z 


L_  /<j,  v  1/2 


¥  r  ©  <«  *  3>-‘ 


rc  ■  Lz(<’i/<’||)1/2  <A1> 


(V'll)' 


2r  /n  +  ML  for.  / <t,, V 
c  zll  llj 


r  »  L 
c 


■z  (V’ll) 


Noting  that  the  critical  radius  is  determined  by  =  T  „  we  find  that  r 

u  cU  c 

is  given  by 


V  L  a,  , 

^2  —  — - -  (M  +  2)  1 

v,  p.  a,, 

l  i  II 


r  <<  L  fa./a..! 
cr  zll  IU 


3-12  Vn  Lz  CT1  .  -.-I 

—  ~  —  77,  (M  *  3> 

11  II 


:r  =  Lz(al/all) 


n  f_n 
2  vv . 


•V  L  cr, 
n  z  1 

,v  .  p .  a. . 

l  i  II 


-z©1/Z) 


r  »  L 
cr  z 


(al/all)' 


where  v.  =  (T  +  T.)/m.  and  p.  =  v,/5?. 

l  e  i  i  l  li 
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It  is  clear  from  (A2)  that  r  is  directly  proportional  to  L  and  V  , 

cr  z  n 

and  inversely  proportional  to  T  ,  1\  ,  and  The  dependence  of  r  on 

1/2 

M  depends  upon  the  value  of  rcr  relative  to  L^ia^/a^)  .  Interestingly, 

1/2 

we  find  that  clouds  are  always  stable  when  r£  >>  L  (o^/<T||) 


and 
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values  are  typical  of  barium  clouds  released  at  -  180  km. 


100  km 
200  km 


values  are  representative  of  a  nuclear  environment. 


1000  m/sec 


1000  m/sec,  670  m/sec,  330  m/sec,  and  100 
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P.o.  BOX  ros 


’EP" 

’ORE.  'LA 

04550 

"K 

..TIN 

DOC 

CON 

FOR 

TECH 

INFO  DEPT 

OR 

ATTN 

DOC 

CON 

FOR 

L-380/R. 

OTT 

O  *  r 

■•TIM 

DOC 

CON 

FOR 

L-  31 

.  R. 

HAGER 

23 


LOS  ALAMOS  NATIONAL  LABORATORY 


i.  BOX 

1663 

;  ALAMOS,  NM 

87545 

01CY 

ATTN 

J.  WOLCOTT 

01CY 

ATTN 

R.F.  TASCHEK 

01CY 

ATTN 

E.  JONES 

01CY 

ATTN 

J.  MALIK 

01CY 

ATTN 

R.  JEFFRIES 

01CY 

ATTN 

J.  ZINN 

01CY 

ATTN 

D.  WESTERVELT 

01CY 

ATTN 

D.  SAPPENFIELD 

LOS  ALAMOS  NATIONAL  LABORATORY 
MS  D438 

LOS  ALAMOS,  NM  87545 


01CY 

ATTN 

S. 

P.  GARY 

01CY 

ATTN 

J. 

BOROVSKY 

SANDIA  LABORATORIES 

P.O.  BOX 

5800 

ALBUQUERQUE,  NM 

87115 

01CY 

ATTN 

W. 

BROUN 

01CY 

ATTN 

A. 

THORNBROUGH 

01CY 

ATTN 

T. 

WRIGHT 

01CY 

ATTN 

D. 

DAHLGREN 

01CY 

ATTN 

3141 

01CY 

ATTN 

SPACE  PROJ  DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

01CY  ATTN  B.  MURPHEY 

01CY  ATTN  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  DC  20545 

01CY  ATTN  DR.  YO  SONG 

NATL.  OCEANIC  &  ATMOSPHERIC 
ADMINISTRATION 
ENVIRONMENTAL  RESEARCH  LABS 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

01CY  ATTN  R.  GRUBB 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 

P.O.  BOX  92957 

LOS  ANGELES,  CA  90009 


01CY 

ATTN 

I. 

GARFIJNKEL 

01CY 

ATTN 

T. 

SALMI 

01CY 

ATTN 

V. 

JOSEPHSON 

01CY 

ATTN 

S. 

BOUER 

01CY 

ATTN 

D. 

OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

01CY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOCIATION,  INC. 
1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

01CY  ATTN  L.  SLOAN 

01CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY.  CA  94701 


01CY 

ATTN 

J. 

WORKMAN 

01CY 

ATTN 

C. 

PRETTIF. 

01CY 

ATTN 

S. 

BRECHT 

BOEING  COMPANY. 
P.O.  BOX  3707 

THE 

SEATTLE, 

VA  981 

24 

01  CY 

ATTN 

G . 

KEISTER 

01CY 

ATTN 

D. 

MURRAY 

01CY 

ATTN 

G. 

HALL 

01CY 

ATTN 

J. 

KENNEY 

CHARLES  STARK  DRAPER  LABORATORY 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  92139 

01CY  ATTN  D.B.  COX 

01CY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
223O0  COMSAT  DRIVE 
CLARKSBURG.  HD  20871 
01 CY  ATTN  G.  HYDE 


CORNELL  UNIVERSITY 

DEPT  >'!■'  EIECTRIcaI.  ENGINEERING 

ITHACA.  NY  1  ..N-n 

"lev  •T~N  i'.E.  FARLEY.  JR. 
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ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

01CY  ATTN  H.  LOGSTON 

01CY  ATTN  SECURITY/ 

(PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  WILSHIRE  BLVD. 

SANTA  MONICA,  CA  90401 

01CY  ATTN  C.G.  GABBARD 

01CY  ATTN  R.  LELEVIER 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

01CY  ATTN  SECURITY  OFFICER 
01CY  ATTN  T.N.  DAVIS 
01CY  ATTN  NEAL  BROUN 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP- 
EASTERN  DIVISION 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
FJDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

01CY  ATTN  DAN  MCCLELLAND 

01CY  ATTN  K.  YEH 


INSTITUTE  FOR  DEFENSE  ANALYSIS 
1801  NO.  BEAUREGARD  STREET 


ALEXANDRIA,  VA 

22311 

01CY 

ATTN 

J.M. 

AEIN 

01CY 

ATTN 

ERNEST  BAUER 

01CY 

ATTN 

HANS 

VOLFARD 

01CY 

ATTN 

JOEL 

BENGSTON 

INTL  TELL  6.  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY.  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 


JAYCOR 

P.O.  BOX  85154 
11011  TORREYANA  ROAD 
SAN  DIEGO,  CA  92138 

01CY  ATTN  N.T.  GLADD 

01CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOC  LIBRARIAN 

01CY  ATTN  THOMAS  POTEMRA 

01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORPORATION 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET 
(P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 


01CY 

ATTN 

DASIAC 

01CY 

ATTN 

WARREN  S.  KNAPP 

OlCY 

ATTN 

WILLIAM  MCNAMARA 

01CY 

ATTN 

B.  GAMBILL 

LINKABIT  CORPORATION 

10453  ROSELLE 

SAN  DIEGO.  CA  92121 

01CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 

01CY  ATTN  DEPT  60-12 

01CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  HISSILES  f.  SPACE  CO.,  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 


OlCY 

ATTN 

MARTIN  WALT/ 

DEPT  52-12 

OlCY 

ATTN 

W.L.  IMHOF/ 

DEPT.  52-12 

OlCY 

ATTN 

RICHARD  G.  JOHNSON/ 
DEPT.  52-12 

ulry 

ATTN 

J  .  B .  C  LA  D I S  / 

DEPT.  52-12 
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MARTIN  MARIETTA  CORPORATION 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01CY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  02647 


01CY 

ATTN 

N.  HARRIS 

01CY 

ATTN 

J.  MOULE 

01CY 

ATTN 

GEORGE  MROZ 

01CY 

ATTN 

W.  OLSON 

01CY 

ATTN 

R.W.  HALPRIN 

01CY 

ATTN 

TECHNICAL  LIBRARY 
SERVICES 

MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  03101 


01CY 

ATTN 

P.  FISCHER 

01CY 

ATTN 

W.F.  CREVIER 

01CY 

ATTN 

STEVEN  L.  GUTSCHE 

01CY 

ATTN 

R.  BOGUSCH 

01CY 

ATTN 

R.  HENDRICK 

01CY 

ATTN 

RALPH  KILB 

01CY 

ATTN 

DAVE  SOULE 

01CY 

ATTN 

F.  FAJEN 

01CY 

ATTN 

M.  SCHEIBE 

01CY 

ATTN 

CONRAD  L.  LONGMIRE 

01CY 

ATTN 

B.  WHITE 

01CY 

ATTN 

R.  STAGAT 

01CY 

ATTN 

D.  KNEPP 

01CY 

ATTN 

C.  RINO 

MISSION  RESEARCH  CORPORATION 
1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NM  87106 

01CY  ATTN  R.  STELLINGWERF 

01CY  ATTN  M.  ALME 

01CY  ATTN  L.  VRIGHT 

MITRE  CORPORATION 
WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

01 CY  ATTN  V.  HALL 

01CY  ATTN  V.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD 
LOS  ANGELES,  CA  90025 

01CY  ATTN  E.C.  FIELD.  JR 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
UNIVERSITY  PARK,  PA  16802 
01  CY  ATTN  IONOSPHERIC 
RESEARCH  LAB 


PHOTOMETRICS,  INC. 
4  ARROW  DRIVE 
WOBURN,  MA  01801 


01CY 

ATTN 

IRVING  L.  KOFSKY 

PHYSICAL 

DYNAMICS,  INC. 

P.O.  BOX 

10367 

OAKLAND, 

CA  04610 

01CY 

ATTN 

A.  THOMSON 

R  &  D  ASSOCIATES 

P.O.  BOX 

9695 

MARINA  DEL  REY, 

CA  90291 

01CY 

ATTN 

FORREST  GILMORE 

01CY 

ATTN 

W.B.  WRIGHT,  JR 

01CY 

ATTN 

W.J.  KARZAS 

01CY 

ATTN 

H.  ORY 

01CY 

ATTN 

C.  MACDONALD 

01CY 

ATTN 

BRIAN  LAMB 

01CY 

ATTN 

MORGAN  GROVER 

RAYTHEON  CORPORATION 
528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42ND  STREET 
NEW  YORK,  NY  10036 

01CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
10260  CAMPUS  POINT  DRIVE 
SAN  DIEGO,  CA  92121-1522 
01CY  ATTN  L.M.  LINSON 

01CY  ATTN  D. A.  HAMLIN 

01CY  ATTN  E.  FRIEMAN 

01CY  ATTN  E.  A .  STRAKEP, 

<3 1CV  ATTN  C .  A .  SMITH 

SdEIirE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
171"  Gu"DP  IL'GE  DRIVE 
HOLE.'.:’.  "A  22102 

"1 1  Y  ATTN  J.  COCKAYNE 

"ICY  ATTN  E.  HYMAN 


26 


SRI  INTERNATIONAL 
333  RAVENSUOOD  AVENUE 


MENLO  PARK,  CA 

94025 

01CY 

ATTN 

J. 

CASPER 

01CY 

ATTN 

DONALD  NEILSON 

01CY 

ATTN 

ALAN  BURNS 

01CY 

ATTN 

G. 

SMITH 

01CY 

ATTN 

R. 

TSUNODA 

01CY 

ATTN 

D. 

A.  JOHNSON 

01CY 

ATTN 

W. 

G.  CHESNUT 

01CY 

ATTN 

C. 

L.  RINO 

01CY 

ATTN 

WALTER  JAYE 

01CY 

ATTN 

J. 

VICKREY 

01CY 

ATTN 

R. 

L.  LEADABRAND 

01CY 

ATTN 

G. 

CARPENTER 

01CY 

ATTN 

G. 

PRICE 

01CY 

ATTN 

R. 

LIVINGSTON 

01CY 

ATTN 

V. 

GONZALES 

01CY 

ATTN 

D. 

MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 


01CY 

ATTN 

U.P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 

ONE  SPACE 

PARK 

REDONDO  BEACH, 

CA  90278 

01CY 

ATTN 

R.K.  PLEBUCH 

01CY 

ATTN 

S.  ALTSCHULER 

01CY 

ATTN 

D.  DEE 

01CY 

ATTN 

D.  STOCKUELL/ 

SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 


BURLINGTON,  MA 

01803 

01CY 

ATTN 

W. 

REIDY 

01CY 

ATTN 

J. 

CARPENTER 

01CY 

ATTN 

C. 

HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01CY  ATTN  N.  ZABUSKY 
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